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We initiate the study of equilibration rates of strongly coupled quark-gluon plasmas in the ab¬ 
sence of conformal symmetry. We primarily consider a supersymmetric mass deformation within 
Af — 2* gauge theory and use holography to compute quasinormal modes of a variety of scalar opera¬ 
tors, as well as the energy-momentum tensor. In each case, the lowest quasinormal frequency, which 
provides an approximate upper bound on the thermalization time, is proportional to temperature, 
up to a pre-factor with only a mild temperature dependence. We find similar behaviour in other 
holographic plasmas, where the model contains an additional scale beyond the temperature. Hence, 
our study suggests that the thermalization time is generically set by the temperature, irrespective 
of any other scales, in strongly coupled gauge theories. 


Introduction and summary.- The discovery of the 
quark-gluon plasma (QGP) in ultrarelativistic heavy ion 
collisions mm is one of the landmark breakthroughs of 
the 21 st century physics. Its behaviour as a strongly- 
coupled liquid made it clear that our ability to describe 
these experiments in an ab initio manner directly from 
quantum chromodynamics (QCD) is insufficient. In par¬ 
ticular, one outstanding theoretical puzzle is elucidating 
microscopic mechanisms behind the rapid appearance of 
the QGP itself (with temperature of order of 500 MeV 
[3]) in less than 1 fm |4|. 

The absence of reliable theoretical tools to describe 
time-dependent phenomena in QCD away from the 
regime of small coupling brought a great deal of attention 
to strongly-coupled gauge theories solvable using holog¬ 
raphy ( e.g., see [5], for a recent review). In holography, 
the QGP creation is associated with the formation of a 
black hole [6], and the QGP equilibration with the relax¬ 
ation of the dual black hole mm. In the latter, equili¬ 
bration rates are associated with imaginary parts of the 
dual black hole quasinormal modes (QNM) frequencies. 

Extensive studies of the simplest holographic mod¬ 
els in which gauge theories are conformally invariant, 
most notably of Af = 4 super Yang-Mills theory (SYM), 
found near- and far-from-equilibrium relaxation times 
not greater than the inverse of the final temperature 
{e.g., see 0 0). As the phenomenological estimates 
obey 1 fm x 500 MeV = 0( 1), one may conclude that 
the large value of the QCD coupling constant is the factor 
responsible for fast equilibration of the QGP. However, 
this conclusion is premature, as, at least near equilib¬ 
rium, the temperature provides the only scale in the A f 
= 4 SYM plasma and dimensional analysis dictates that 
the equilibration rate is proportional to the temperature. 

Motivated by this crucial point, our letter considers 
equilibration processes in the holographic model describ¬ 
ing the nonconformal Af = 2* gauge theory. This theory 
is obtained as a (supersymmetric) deformation of Af = 4 
SYM, giving mass m to half of its fields [TO] . The QGP 
plasma of Af = 2* gauge theory is characterized by the 


temperature T and one dimensionless parameter m/T , 
the ratio of the mass to the temperature. We vary the 
latter from m/T = 0 {i.e., the Af = 4 SYM plasma) up 
to m/T = 25. The biggest deviation from conformality 
occurs at m/T ~ 4.8 (see Fig. 0 m • 

We primarily focus on the behaviour of small homo¬ 
geneous perturbations of AT = 2* plasma given by phe¬ 
nomenological operators of dimensions 2, 3 and 4, top- 
down operators of dimensions 2 and 3 and the energy- 
momentum tensor m We compute the frequencies of 
the corresponding least-dampened QNM, as they set an 
approximate upper bound on thermalization times. We 
observe that QNM frequencies, when normalized to tem¬ 
perature, change rather mildly as a function of m/T (see 
Fig. |2j [4] and [ 5 ]) . The largest change in frequencies with 
respect to Af = 4 SYM is roughly a factor of two and is 
typically much less. This implies that the equilibration 
time for generic small perturbations is still set by T, as 
in the conformal case 0. 

To test the robustness of this finding, we also consid¬ 
ered other nonconformal holographic models. While we 
do not present the details, examining the (nonsupersym- 
metric) bosonic mass deformation of Af = 4 SYM [10] 
yields qualitatively similar behaviour for all of the above 
quantities. We also consider QNM of a strongly coupled 
conformal plasma with chemical potential p. This con¬ 
stitutes another example of equilibrium state specified by 
a single dimensionless number, i.e., p/T. We analyze the 
behaviour of phenomenological operators of dimensions 
2, 3 and 4. We find that there is a regime in which the 
lowest QNM frequencies, when normalized to tempera¬ 
ture, change mildly as a function of p/T (see Fig. |6|. 

Finally, let us notice that there is strong evidence that 
equilibration in 1/T also occurs for large deviations from 
equilibrium mHm. Here, our results in Fig.[3]show that 
including momentum dependence does not change the 
qualitative behaviour of Af = 2* QNM, as compared to 
Af = 4 SYM. Moreover, the so-called close limit approxi¬ 
mation mm, originally devised for black hole mergers 
in asymptotically flat spacetimes, suggests that nonlinear 
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effects as seen by the expectation values of local opera¬ 
tors in a boundary field theory are not decisive once a 
single horizon forms in the bulk PE]. We are thus led 
to conjecture that equilibration in any strongly coupled 
gauge theory, once a single horizon forms in a dual grav¬ 
ity description, will be captured in a satisfactory way 
by Af = 4 SYM (see [2(Tl [21] for evolution of bottom- 
up nonconformal holographic models before this occurs). 
For holographic descriptions of heavy ion collisions, this 
suggests the most fruitful direction will be relaxing geo¬ 
metric symmetry constraints in the dynamics in Af = 4 
SYM (e.g., see [22 HU) rather than writing complicated 
codes to incorporate broken conformal symmetry. 

Af = 2* thermodynamics.- The holographic descrip¬ 
tion of Af = 2* QGP is a five-dimensional black hole 

ds 2 = e 2A M (-B(z) 2 dt 2 + dx 2 ) +L 2 ^, (1) 

where L sets the overall curvature scale and the boundary 
is at z = 0. The metric 0 solves the equations of motion 
derived from the following effective action 

s = J d 5 ^ %/-£ (R + ^matter) , (2) 

where the matter part, as discovered by Pilch and Warner 
in [24] (PW), is given by 

- ^matter = —£>pw 12(<9<a) 2 + 4(<9y) 2 + V (3) 

and the scalar potential takes the form 

V = T [~4e~ 4a - 8e 2a cosh 2 X + e 8a sinh 2 2 X ] ■ (4) 

±J Z 

The mass m in the Af = 2* theory is encoded in the 
near-boundary behaviour of the scalars as 

X = ^z + Oiz 2 log z), a= kT z 2 log z + 0(z 2 ). 

( 5 ) 

We choose the units in which the event horizon is located 
at z = 1 (i.e., B(l) = 0). The QGP temperature is then 
given by T = e A ^ B' (1)/ (2ttL) and there exists a one- 
parameter family of geometries labeled by the ratio m/T. 
We construct those geometries using a shooting method. 
See mmm f° r details of Af = 2* holography, its equi¬ 
librium thermodynamics and numerical implementation 
of the shooting method. 

As a benchmark of Af = 2* thermodynamic properties, 
in Fig. [l]we consider the trace of the stress tensor nor¬ 
malized to the energy density of Af = 4 SYM and plot 
this ratio as a function of m/T. This quantity starts at 
0 for m/T = 0, as Af = 4 SYM is conformally invariant, 
and peaks at m/T « 4.8. We thus expect the strongest 
violations of conformality to occur precisely there. 

Thermalization time for phenomenological scalar 
operators.- The simplest equilibration process to con¬ 
sider is perturbing a holographic QGP by exciting a phe¬ 
nomenological operator O of dimension A. The gravity 
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FIG. 1. Trace of the energy-momentum tensor normalized to 
the energy density of Af = 4 SYM (eo = §7 t 2 N/T 4 with N c 
denoting the number of colours) as a function of m/T. The 
results indicate that, thermodynamically, the effects of the 
conformal symmetry breaking are the strongest at m/T w 4.8. 
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FIG. 2. Real (green continuous) and minus imaginary (red 
dashed) parts of the lowest quasinormal mode frequencies for 
phenomenological scalar operators of dimensions A = 2, 3 and 
4 (from bottom to top). The frequencies do not change sig¬ 
nificantly as a function of m/T, which leads to universal equi¬ 
libration in 1/T. One can also infer from this plot that all the 
frequencies asymptote at low temperatures to the quasinormal 
mode of a massless scalar field living in the (l+5)-dimensional 
AdS-Schwarzschild geometry (dotted curves). For explana¬ 
tion, see the supplemental material. 


description of the latter is a minimally coupled scalar if: 

s = J d 5x V^g{{dip) 2 - m2 iP 2 ) ( 6 ) 

where the mass is given by M 2 = A(A — 4)/L 2 . 

We search for solutions with a harmonic dependence 
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FIG. 3. Momentum dependence of the real (green) and minus 
imaginary part (red) of the QNM frequency of phenomeno¬ 
logical operators with A = 2, 3 and 4 (from bottom to top) 
for m/T = 0 (A f — 4 SYM, dashed) and m/T = 4.8 (con¬ 
tinuous). Surprisingly, corresponding curves are very close to 
each other despite of the fact that m/T = 4.8 matches the 
locus of the maximal deviation from conformal invariance in 
thermodynamics of Af = 2* (see Fig. [I]). 


general, monotonic functions of m/T. For an explana¬ 
tion of why all the curves in Fig. [2] are saturate in the 
low temperature regime, see the supplemental material. 

In Fig. [3] we set m/T = 4.8 where (approximately) the 
biggest deviations from conformality are expected and 
plot A = 2,3 and 4 QNM frequencies as a function of 
momentum normalized to temperature. We find not only 
that deviations from the Af = 4 SYM result are small, but 
also that up to an overall normalization, the momentum 
dependence is almost the same in all three cases. 

Thermalization time for physical Af = 2 * gauge theory 
scalar operators.- We would like to extend analysis of 
the previous section to physical operators of the Af = 2 * 
gauge theory. In [29], Bobev, Elvang, Freedman and Pufu 
(BEFP) showed that agreement between gravity calcula¬ 
tions and exact QFT results in the J\f = 2* gauge theory, 
obtained using localization techniques m , requires em¬ 
bedding of the PW action © into a larger consistent 
truncation of five-dimensional Af = 8 gauged supergrav¬ 
ity. The matter part of the action © then becomes 


(dp) 2 

— £ matter = ~£bEFP = 12-x-h 


4(<9X ) 2 
(1 -X 2 ) 2 


+ V (7) 


where the potential V is 


on t and x, i.e., if = e~' luJt+lkx f {z). We use the shooting 
method to search for the normalizable modes satisfying 
an ingoing boundary condition at the horizon. These 
conditions are satisfied by a discrete set of complex fre¬ 
quencies that define QNM. In Fig. [2] we set k = 0 and 
plot the real and imaginary parts of the lowest QNM 
frequencies normalized by 2eT as a function of param¬ 
eter m/T for A = 2,3 and 4. One can clearly see that 
the frequencies do not change significantly {i.e., by more 
than a factor of 2) as a function of m/T. This results 
in universal thermalization time set by the inverse of the 
temperature. Note that the QNM frequencies are not, in 


l + Vi±X_ rfi Tl ) 2 + (*») 2 ' 

if l-X 2 (1-V 2 ) 2 

(8) 

The vector X = (Xi, X 2 , X 3 , X 4 , X 5 ) combines five of the 
scalars and 77 is the sixth. The symmetry of the action 
reflects the symmetries of the dual gauge theory [29] : the 
two scalars (Xi,X 2 ) form a doublet under the U{1)r, 
while (X 3 , X 4 , X 5 ) form a triplet under SU{2)y and 77 is 
a singlet. The PW action © from 0 with 

e a = p, X\ — tanh x ? X 2 = X 3 = X 4 = X 5 = 0 . (9) 



The effective action describing fluctuations of PW backgrounds within BEFP is obtained linearizing 0 in 
X 2 ,X 3 ,X 4 ,X 5 scalar fields: 


SC = —Cbefp + £pw + = SC 2 + SCy , 

SC 2 = -(1 + c) 2 (dx 2 f - (V + c)pg /3 - 4(1 + c)pl /3 + ff/ff) (X2) 2 , 

5Cy = -(1 + cf(dXy ) 2 - FF { {c 2 _ 1)/9 4/3 _ ^ + ^1/3 + ) {jty) 2 (10) 


where Xy = (X 3 ,X 4 ,Xs) and p§ = e 6a , c = cosh2%. 
Since SC is SU{2)y invariant, it is enough to consider 
only one component of Xy, say X 3 . The scalars X 2 and 
X 3 correspond to operators 0 3 and 0 2 with dimensions 


A = 3 and 2 of Af = 2 * gauge theory. Following [31], we 
compute the lowest QNM frequencies of these operators. 
The results in Fig. [4] show once again that the QNM 
frequencies have a very mild dependence on m/T. 
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FIG. 4. Real (green continuous) and minus imaginary (red 
dashed) parts of the lowest quasinormal mode frequencies for 
O 3 and O 2 operators (see (10)) in J\f — 2* gauge theory (from 
top to bottom). The frequencies do not change significantly 
as a function of m/T. 
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FIG. 5. Real (green continuous) and minus imaginary (red 
dashed) parts of the lowest quasinormal mode frequencies the 
operators O x and O a (from top to bottom) in AT = 2* gauge 
theory. 


The matter fluctuations 5a and 5x represent operators 
with A = 2 and 3 in the J\f = 2* gauge theory. The 
results for the lowest QNM frequencies of 5a and 5\ are 
collected in Fig. [5] As in the previous examples, the 
frequencies exhibit a very mild dependence on m/T. 

Thermalization time for the energy-momentum 
tensor.- Quasinormal modes of the energy-momentum 
tensor in the boundary theory are dual to metric pertur¬ 
bations in the bulk. Different polarizations decompose 
into decoupled sets of various helicities with respect 
to the propagation direction in the boundary [9]. The 


simplest helicity-2 fluctuations are always equivalent 
to a minimally coupled massless scalar [32] and its 
momentum-dependence (at m/T = 4.8) is given by 
A = 4 curve in Fig. i ln the zero-momentum limit, 
the helicity distinction vanishes implying that generic 
metric perturbations obey the equation of motion for a 
minimally coupled massless scalar. Hence in this case, 
the QNM are the same as the scalar operator with 
A = 4 plotted in Fig. [2] Finally, in Fig. [5] we plot 
zero-momentum QNM associated with the perturbations 
of the background scalars a and Again, we find rather 
mild dependence on m/T. 

Comparison to conformal plasma at nonzero chemical 
potential.- From an operational point of view, deviations 
from conformality in J\f = 2* gauge theory amount to the 
presence of a dimensionless parameter m/T characteriz¬ 
ing the equilibrium configuration. In order to confirm the 
robustness of our main result, we consider a strongly cou¬ 
pled conformal QGP at nonzero chemical potential fi. In 
this case, the dimensionless parameter is pb/T. The dual 
solution is the anti-de Sitter Reissner-Nordstrom black 
hole (AdS-RN). We compute the lowest zero-momentum 
quasinormal modes of phenomenological operators (with 
A = 2, 3 and 4) dual to probe (neutral) scalar fields in 
the AdS-RN background. 

Again, as shown in Fig. [6] the frequencies depend 
mildly on (i/T unless the chemical potential (at fixed 
temperature) gets big enough to see the effects of the 
critical behaviour associated in the dual description with 
the (near-)extremal AdS-RN. In this case, the dominant 
QNM has a purely imaginary frequency. We explain in 
the supplemental material why we do not expect similar 
transition to occur in M — 2* gauge theory and refer the 
reader to refs [33: and [34], which discuss these purely 
imaginary QNM in detail. 

Note added: While this letter was being finalized we 
learnt about the upcoming results of Ref. [35], which 
presents the computation of QNM in a bottom-up holo¬ 
graphic QCD model. The results of f35] are in line with 
the point of view presented in this letter, i.e. equilibra¬ 
tion rates in strongly coupled nonconformal plasma are 
not very different from those of M = 4 SYM. 
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draft. We would like to thank the authors of articles m 
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cation. This work was partly supported by the National 
Science Centre grant 2012/07/B/ST2/03794. Research 
at Perimeter Institute is supported by the Government 
of Canada through Industry Canada and by the Province 
of Ontario through the Ministry of Research & Innova¬ 
tion. AB and RCM acknowledge support from NSERC 
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FIG. 6. Real (green continuous) and minus imaginary (red 
dashed) parts of the QNM frequencies originating at p/T = 0 
from the lowest QNM for phenomenological scalar operators 
of dimensions A = 2,3 and 4 (from bottom to top) in the 
AdS-RN background. Red dot-dashed curves denote purely 
imaginary QNM for operators of dimensions A = 2, 3 and 4 
(also from bottom to top), which are special to finite p and for 
sufficiently large fi/T become the least-dampened modes (see 
(33J [34] for the discussion of this branch). Before this transi¬ 
tion occurs, the QNM frequencies change by a small amount, 
similar to what we see for AT = 2* gauge theory. Despite the 
transition, the equilibration time remains to be bounded by 
1/T, in accord with the main message of this letter. Finally, 
note that it is not possible to relate in any meaningful way the 
values of parameter fi/T appearing here and m/T of JV = 2* 
gauge theory, as they describe two different physical systems. 
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FIG. 7. Trace of the energy-momentum tensor normalized to 
the energy density as a function of m/T (see also Fig. [TJ . This 
quantity saturates at low temperatures to the value consistent 
with the thermodynamics of a five-dimensional CFT. 


red (IR) region of the five-dimensional geometry 
2 / 9 \ ^ 3 

e2 *“ 2 “' ^ ~ ’ 'MS?) (mL> • (11) 

where we introduced a new radial coordinate u. Defining 
4>X — ^ (3a + x) and </>2 = 1 (a - x) (12) 
one can lift the bulk geometry 0 to six dimensions with 
dsl = e~ 2 ^ 2 ds\ + 4 m 2 L 2 e°^ 2 dx\ . (13) 


Supplemental material- Ref. [36], following 
observed that for large m/T the thermodynamics 
and the lowest transport coefficients of the four¬ 
dimensional AT = 2 * gauge theory are inherited from a 
five-dimensional conformal field theory (CFT) compacti- 
fied on a circle. This is apparent in Fig. [7| where the ratio 
of the trace of the thermal stress tensor to the energy den¬ 
sity approaches 0.25 for large m/T. In a five-dimensional 
CFT, the energy density and pressure are related by 
65 = 4P 5 , which indeed gives (es — 3 P$) /c^ = 1/4. 

Below we will demonstrate that this near-horizon AdS 
region dominates also the behaviour of the QNM in the 
regime of large m/T. This constitutes a highly nontriv¬ 
ial test for our numerics and demonstrates that we did 
not miss any exotic branches of QNM’s. To set the back¬ 
ground, let us thus briefly remark on the relation between 
the five-dimensional bulk geometry in the low tempera¬ 
ture regime ( m/T 1) and an auxiliary six-dimensional 
dual of a finite-temperature CFT [36]. At low enough 
temperature the field x attains a large value in the infra- 


In particular, the zero-temperature geometry recovered 


in the IR using (13) is a six-dimensional AdS spacetime 
3 3 / 2 !/ 2 


ds 51 ~ 


2 u 2 


AfYl ^ 

du 2 H—— ( rin^dx^dx v + dxf) 


(14) 


where the six-dimensional curvature scale is Lq/L = 
3 3 / 4 /2 1 / 2 . 

Let us now reconsider the QNM of phenomenological 
operators from the point of view of the six-dimensional 
geometry in the IR. Uplifting the scalar field action 0 
to six dimensions using (13) gives 


S ~ j d 6 xx^—§ ((dip) 2 + e 2 ^ 2 M 2 ip 2 ^j , (15) 

where tilde denotes quantities defined using the six¬ 


dimensional metric (13). As in the relevant regime 


approaches — oo, see eqs. ( 11 ) and ( 12 ), the mass term 


in eq. (15) becomes negligible and the resulting QNM 


frequencies must be that of a massless scalar field in a 
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six-dimensional AdS-Schwarzschild black hole. Recalling 
the holographic dictionary for scalars in AdS6/CFT 5 , 

M 2 = ^ A(A-5), (16) 

we see the massless scalar in the effective six-dimensional 
geometry is dual to a marginal operator in the dual five¬ 
dimensional CFT. This is in excellent agreement with 
our numerics, as fits to the large-ra/T behaviour of fre¬ 
quencies plotted in Fig. [2] are all consistent with the 
lowest QNM of a marginal operator in a five-dimensional 
holographic CFT. We performed similar analysis for the 
other fields which we considered and found analogous be¬ 
haviour with the exception that the effective IR masses in 
uplifted actions need not always be zero. Let us also note 
that in Fig. [4] and [5] the frequencies saturate to their IR 
values much more slowly than in Fig. [2] because of com¬ 
plicated couplings of the corresponding perturbations to 
the background fields. 

This result is interesting on its own right. To the best 
of our knowledge, it provides the first example in hologra¬ 
phy where the low temperature non-equilibrium response 
(i.e.j away from uj -» 0) of strongly coupled matter is fully 
determined by the properties of the IR geometry in the 
dual gravity description. 
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